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Abstract 

A result from Dodd and Gibbs]l] for the second virial coefficient of particles in 1 dimension, 
subject to delta-function interactions, has been obtained by direct integration of the wave func- 
tions. It is shown that this result can be obtained from a phase shift formalism, if one also includes 
the contribution of oscillating terms. The result is important in work to follow, for the third virial 
coefficient, for which a similar formalism is being developed. We examine a number of fine points 
in the quantum mechanical formalisms. 

Introduction 

The quantum mechanical formulation of virial coefficients, or the calculation of corrections to classical 
results, continues to be of concern. We note, for example, the fairly recent article by Hussien and Yahia 
P], on the quantum corrections to the classical fourth virial coefficient for square well potentials. Our 
efforts are, instead, directed to the continuing development of a fully quantum mechanical formalism, 
applicable to the higher virials, and especially useful at low temperatures. 

We recall that, many years ago, Dodd and Gibbs^lJ were able to evaluate the second and third 
virial of identical particles in one dimension, subject to repulsive delta function interactions. To do 
this, they introduced a complete set of energy eigenfunctions in the traces, and then integrated to 
obtain answers in terms of single and two-dimensional integrals. 

Their original aim, as is clear from their paper and from Gibbs' thesis, was to obtain a formalism 
in terms of the S-matrix, following the formalism of Dashen, Ma and Bernsteinpl, in an effort to 
obtain the virials in terms of scattering quantities. This was not possible, as the authors determined 
that the 'formal limiting processes are not valid for the singular amplitudes of this system'. 

An alternative formalism, proposed by Larsen and MascheroniQ, aimed at expressing the higher 
virials in terms of eigenphase shifts and bound state energies, and, if necessary, other on-shell prop- 
erties, in a hyperspherical harmonic formalism. Under some constraints, one is then able to recover 
the classical results for the 3rd virial, through a semi-classical approximation^- Low temperature 
calculations have been carried out for repulsive potentials in 2 dimensions and a more sophisticated 
formalism, involving a hyperspherical adiabatic basis, has been presented|7], 

Our aim, in this and subsequent papers, is to use the excellent results of Dodd and Gibbs to confirm 
our virial formalism for the second and third virials, for this particular dimension and interaction, for 
which we have analytical and numerical resultsjBJ El HOI- More specifically, in this short note, we will 
show that, for the second virial, the result of Dodd and Gibbs can be rewritten as an integral, over 
the energy, of a Boltzmann factor together with the phase shift - the usual second virial type of result 
- plus that of an extra 'oscillating' type of term, particular to the contribution of even parity wave 
functions in one dimension. This situation is understood^!], though not necessarily widely, and we 
discuss it for this particular case and in greater generality. 



1 



The Second Virial Coefficient 



Following ref. [Jj, in this section we present the main steps leading us from the grand partition function 
to the expression of the second virial coefficient in terms of the asymptotic properties of the system. 
Instead of using the Ursell expansion, we simply take, directly, the logarithm of the grand partition 
function. 

The latter can be written as 

Q = j>2 z n Tr ( e -^») = 1 + zTr (e^ 1 ) + z 2 Tr (e^ 2 ) + • • • (1) 
n=0 

The fugacity z equals exp (jti/kT); (3 = 1/kT, where [i, k and T are the Gibbs' function per particle, 
Boltzmann's constant and the temperature, respectively. H n and T„ are the n-particle Hamiltonian 
and kinetic energy operators. 
It is important to note that 

Tr (e _/3ff ") ~>V n asV -> large. 
However, the logarithm of the partition function 



InQ = ln(l + \zTr(e-P Tl ) + z 2 Tr{e~^ H A + ■■■]) 
[zTr (e~^ )+z 2 Tr (e"^ ) + •••] 
-\ [zTr (e-^) + z 2 Tr (e"^) + • ■ ■] 2 

or, equivalently, 
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(2) 



when divided by V, gives coefficients which are independent of the volume, when the latter becomes 
large; we call them b\. We can then write for the pressure and the density 



V V z . 



where 



h = ^Tr (e-^) ; b 2 = i {rr (e^) - \ [Tr (e"^)] 2 } . etc. 
Elimination of z from Eqs. © gives the virial expansion of the equation of state 



kT ^ ' 



P ■ 



Let us evaluate the fugacity coefficients. 

bl =±Y / < k \e~ {h2/2m)0k2 \ k > 
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V (2^F 

1 _i_ 

(h 2 /2irmkT) 3 / 2 ~ 



(3) 



(4) 



(5) 



2 



where we introduced a complete basis set, for positive energies in the continuum, and wrote the result 
in terms of the thermal wave length, At- 

Now, let us to evaluate the second coefficient in three dimensions. Consider the Hamiltonian H 2 , 
separated in a CM part and a relative one, H 2 = H^" 1 + H% . Note that for the center of mass 



to -> 2m, thus A T 3 -> 2 3/2 A T 3 . 

The trace of the CM part becomes 



^ r (e-^~)=Ir,( e -«-) = ^ (6) 



and the trace of the relative part, in the symmetric case, can be written as 

2V (V^ re! ) = < m\e-^ e '\m > + J dq < q\e~^\q > 

m 
m 

where we assumed a spherically symmetric potential and introduced a complete set of symmetric 
eigenstates. The index q includes the continuum variable k, as well as the discrete angular momentum 
indices I and to. The integration over "r* is carried out over the reduced range (half) of u?, since 
this then yields all of the configurations of the symmetrical wave functions. These wave functions are 
normalized over this reduced range. Our result reads: 



poo />oo 

Yfe-f > <+ / dke-P E >< dr 

™ l = even *'° 



Xkt(r)Xkl(r) ■ (8) 



A method, that has been used to evaluate the integral over k, is to put the system in a box of 
volume V (see, for instance, ref. 0j). The energies will then be discrete, characterized by an index n, 
and the eigenfunctions will be normalized so that their integral over the r gives 1. For large volumes, 
we can again proceed from the sum to an integral in which, to estimate the density of states 4|, we 
use the asymptotic form of the wave function 

XM (r) r largc ~ sin ^fcr - y + S e (k) 

We require that it go to zero at the boundary of a large sphere, of radius R, and that therefore the 
argument equal nir. It then follows that 

dn 1 / _ d . . , . 

and that the continuum contribution equals 

Tr s (e~^" - e-^ ! ) = £ (2£ + 1) i J°° dk (±6 e (/.•) ) , 

i — even 

The difference 6| — ^2° m the second fugacity coefficients is then 



e-^ k2 (9) 



2 3/2 

w - 

1 t — even 

and the virial coefficient takes the well known form ^3 m terms of the energies of the bound 
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states and the derivative of the eigenphase shifts: 
B 2 (T) - B\ (T) = 



-N2 3 / 2 A^ 



-PK 



(10) 



An integration by parts can change the integral to involve the phase shifts directly. 



Oscillatory Terms 

The passage from the sum, in the box approximation, to the integral, in the continuum situation, 
is usually perfectly justified. The correction terms are negligible, except under some very specific 
circumstances, as we shall see. 

To investigate this, we follow a method [12] devised for the second virial coefficient for anisotropic 
interactions, where it becomes impossible, in a box, to simultaneously satisfy a common boundary 
condition for all the amplitudes (belonging to different orbital angular momenta) of any given eigen- 
function. It became useful then, and elegant, to directly evaluate the integral of the square of the 
wave functions that already appears (jHJ in our formulation for the symmetric case. To do so requires 
a simple 'Green's function trick'. 

The Schrodinger equation for the radial wave function corresponding to fc is 

" "SrXiW (r) + ^±^X« (r) + ^~V (r) X U (r) = fc 2 Xki (r) . (11) 

Multiplying Eq. ifTTl) by x%'£ ( r )' an d the complex conjugate of the corresponding equation for k' by 
Xu (r), and then taking the difference between them, we obtain 



Xki (r) -rXh'l ( r ) ~ Xk'i 0) jXm (r) 
ar ar 



X*k>i(r)Xkd r ) ■ (12) 



(fc' 2 - fc 2 ) dr 

Integrating over the interval < r < r max , using the boundary conditions at the origin Xki (0) = 0, 
Xk't (0) = and L'Hospital rule, expression l(T2|l becomes 

Yk { {§k Xkl (r) ) t Xh (r) " Xkl (r) Md~r XM (r) L r _ = I drXM {r) XM {r) (13) 
when fc' — > fc. For r max large enough, we can use the asymptotic form of the wave function 

1 



sm(kr-£ir/2 + S e (k)) (14) 

to evaluate the left hand side of the Eq. (flSjl as 

£ (r,nax + -^St (fc)) - ^ sin (kr max - £n/2 + 5 t (fc)) cos (kr max - £n/2 + 5 e (fc)) 

= \ {r max + &8 t (fc) - ± sin [2 (kr max - £ir/2 + S e (fc))] } . [ > 

If we subtract from 11511 the corresponding expression that would be obtained if there were no inter- 
actions, we obtain for the oscillatory terms 

(-1)* 

— 2~Kk~ ^ krmax ^ ( cos t 2 ^ £ (*)] - 1) + cos (2kr max ) sin [25 1 (fc)]} 

which has to be integrated over fc, with a Boltzmann factor. Consequently, what we wish to estimate 
are integrals such as 
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h = JT dk sin[2kr max ] ^iMzT j e -&* 

and 

h = £> dk cos (2kr max )^We-^\ 

In the vicinity of zero Si (k) — (n + e/2)ir + O (k 2l+1 }, where n is the number of bound states 
supported by the potential and e is zero except when I = and there is an s-wave zero-energy 
resonance, in which case e equals 1. Evaluating the integrals, either directly, or making the change 
of variable u = 2kr max , and letting r max become large, we then find that there are no contributions 
to the virial from the oscillatory terms in our three dimensional case, unless there is a zero-energy 
resonance at I — 0, in which case the phase shift starts at ir/2 and I\ equals -it. The total contribution 
for I — 0, to be added to 1/ir times the derivative of the phase shift, is then 1/2. 



The One Dimensional Case 

In this section we discuss the 2nd-virial of a Bose one-dimensional system of particles with repulsive 
^-function interactions. 

The symmetric wave function for the 2-body system J^j m Cartesian coordinates is 

r klM {x u x 2 ) = j= (£) -jM== [(i _ m) e i( fcl * 1+ fc 2 * 2) + (i + m) e i(* 2 * 1+fcl * 2)] (16) 

where Xi is the position coordinate of particle i and k\ and ki are the two momenta (wave numbers) 
involved in the system. The expression is only valid for x 2 > X\, but is correct for all momenta. The 
normalization is such that integrating over these values of Xi gives the expected delta function nor- 
malization, as is appropriate for a symmetric function. The strength of the interaction is represented 
by g; we define g = (mg/2?i 2 ). 

We note that an expression ^H] valid for xi < xi, requires one to interchange the coefficient of the 
exponentials in iflfij) . 

To derive an expression for the second virial coefficient of the one dimensional system, we parallel 
the discussion for the three dimensional case in the earlier sections. 

We begin the discussion of the virial by rewriting (|16(l . and the wave functions in the second of 
Eqs. 0, in terms of center of mass and relative coordinates. We then evaluate the CM. part as in 
Eq. © which, for a system contained on a line of length A, equals 



i Tr s f e -ftr m \ = i r°° dKe -{n-/ im )pK- 

A \ ) 2tt J — oq 

= 2 1 /2/ At , 



(17) 



where K — k± + fe, and the associated position coordinate X — (xi + X2)/2. 
The wave function for the relative motion is then 



r k {x) = -=(-) 



1 , 1 M/2 ^[ 



2! 2tt' y/k 2 + g 
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l _ t JL\ e ikx + 1+ 19\ -ite 



k I \ k 



(18) 



where k — (k 2 — hi) /2 and the relative position is x — x 2 — x\ > 0. 

A representation of this equation, valid for all values of x, and all values of k, can then be written 

as 

r k (x) = (-) 1/2 ™s(k\x\+6(k)), (19) 

7T 

in terms of the phase shift 

5(k) = -- + arctan~, for k > 0, and S(-\k\) = -S(\k\). (20) 
2 g 
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We note that our expression for the wave function is completely symmetric in all of our variables! 
We have the choice of either integrating over half of the x space and the complete k space, or over 
the complete x space and half of the k space. We choose to do the latter. For k positive, we first 
integrate, over x > 0, the square of the wave function. We can do this by using our Green's function 
'trick' as we have done previously for 3 dimensions, but since we know the explicit form of the wave 
function in terms of its phase shift, it is tempting to evaluate the integral directly! 



f 

Jo 



cos 2 (kx + 8{k))dx = —x ri 



1 ■ sin(2<5(fc)) + -!- sm{2kx max + 26{k)), (21) 



/0 2 4fc v v " 4fc 

where we delay letting the upper limit x max go to infinity, until the subtraction of the integral of the 
square of the 'free' solution, and the subsequent evaluation of the 'oscillating' terms, when integrating 
over k. We don't immediately see the term in the derivative of the phase shift, that was obtained 
from the limiting procedure in our former treatment, but, we note from Eq. I|2f)|) that 



Integrating over 
we obtain 



to x T , 



(a simple factor of 2, due to symmetry), and with the normalization, 



LA. 

■k dk K 



(22) 



: S ( k ) + 2^k [sin(2/cx maa; )(cos(2(5(A:)) - 1) + cos(2kx max ) sin(2<5(fc))] . 
We recognize the form of the oscillating terms, which apart from a change of sign, is precisely that 
encountered in three dimensions. 

Following our previous comments on the domain of the integrations, we integrate, over k, from to 
oo, and evaluate h and h, using the phase shift given by (pHf, This gives h = — ir/2, I 2 = 0. Finally, 
joining the contributions of the oscillatory and non-oscillatory terms, we achieve the expression we 
were looking for, for the difference of the fugacity coefficients in terms of asymptotic quantities, 



b° 2 = 



2 l/2 



' dke ^/ m )f^ld5 {k) 
7r dk 



(23) 



Equivalence of the Results of Dodd and Gibbs. 

In this section we show that our formalism, and (f23|) . gives us an equivalent result for the difference 
in the fugacity coefficients, to that obtained in ref. Q]. While the former was derived directly from 
the partition function and is expressed in terms of asymptotic quantities, the latter was derived using 
the usual cluster expansion 22 an d is expressed in terms of integrals. 
The derivative of the phase shift reads 

d5 g 
dk k 2 + g 2 

and, therefore, the integral of the phase shift term in Eq. 1231) is 

(^^I-dke-i^)^^ 
= 2^ ( i )e ( ft *M^ (Verf 



+9 2 



(24) 



(25) 



Joining the contributions of the oscillatory and non-oscillatory terms, we obtain 



1 



1 



= (n 2 /™)/3s 2 



1-erf 



g^](h 2 /m)p 



4 2 2 1 / 2 \ T 
This is exactly the result from Eq.(13) of ref. Q], 

b2 - b 2 = 2ir(r//2 m) Io°° ds eX P [" S V (i h2 / m ) P)] I eX P (~ 2 3 S ) 

when integrated over the variable s . 



(26) 
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